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Reduction  of  air  pollutants  formation  from  hydrocarbon  combustion  process  requires  improvements  in 
combustion  systems.  The  moderate  and  intense  low  oxygen  dilution  (MILD)  combustion  technique  is 
an  opportunity  to  achieve  such  a  goal.  MILD  combustion  is  a  combustion  regime  which  can  be  attained 
by  high  temperature  preheating  and  high  level  dilution.  In  this  paper,  the  mechanism  of  CO  and  C02  for¬ 
mation  for  a  CH4/H2  fuel  mixture  is  studied  under  MILD  combustion  condition  of  a  jet  in  hot  coflow  (JHC) 
burner.  This  investigation  is  done  using  the  computational  fluid  dynamics  (CFD)  and  also  zero  dimen¬ 
sional  well-stirred  reactor  (WSR)  analysis.  The  RANS  equations  with  modified  k-s  equations  are  solved 
in  an  axisymmetric  2D  computational  domain.  The  DRM-22  reduced  mechanism  is  considered  to  repre¬ 
sent  the  chemical  reactions.  The  effects  of  oxidizer  oxygen  concentration  and  fuel  hydrogen  content  are 
studied  on  methane  oxidation  pathways.  Results  show  that  the  higher  hydrocarbon  oxidation  pathways 
are  effective  on  CO  and  C02  formation  under  MILD  condition.  In  the  methane  oxidation  mechanism,  the 
ratio  between  the  main  route  and  ethane  route  is  the  main  reason  of  CO  increment  at  higher  02  level 
under  MILD  condition  in  JHC  laboratory  burner.  The  WSR  analysis  illustrates  that  a  decrease  of  02  con¬ 
centration  in  oxidizer  does  not  necessarily  lead  to  lower  production  of  CO  and  C02. 

©  2013  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Greenhouse  effects  are  results  of  contribution  of  the  burning  of 
fossil  fuels.  Carbon  dioxide  (C02)  is  one  of  the  primary  greenhouse 
gases  which  comes  from  combustion  of  carbonaceous  fuels. 
Although  carbon  monoxide  (CO)  could  contribute  to  greenhouse 
effects  and  is  known  as  an  air  pollutant,  it  is  usually  omitted  when 
discussing  greenhouse  gases  as  its  life  time  is  short  compared  to 
C02.  Focusing  of  emission  standards  on  regulating  pollutants  re¬ 
leased  by  combustion  processes  is  forcing  the  combustion  systems 
toward  more  efficient  technologies.  In  this  framework,  MILD  com¬ 
bustion  is  known  as  a  promising  technique  to  eliminate  the  rate  of 
pollutant  generation  in  combustion  systems.  Higher  thermal  effi¬ 
ciency  with  this  technique  leads  to  lower  fuel  consumption  which 
demonstrates  suppression  of  C02  generation.  Furthermore  a  more 
complete  combustion  process  limits  the  CO  contribution.  MILD 
combustion  is  acronym  of  ‘Moderate  or  Intense  Low-oxygen  Dilu¬ 
tion’  combustion.  MILD  combustion,  HiTAC,  HiCOT,  and  flameless 
combustion  have  been  broadly  known  as  same  technologies.  This 
combustion  regime  has  been  achieved  using  both  high  tempera¬ 
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ture  preheating  of  reactants  and  high  degree  dilution  of  reaction 
zone,  in  such  a  way  the  temperature  increase  in  the  combustion 
process  is  lower  than  the  mixture  auto  ignition  temperature  [1]. 
Several  numerical  and  experimental  investigations  on  this  com¬ 
bustion  mode  have  shown  specific  characteristics  of  this  new  com¬ 
bustion  regime.  These  characteristics  are  briefly  lower  emission 
production  and  fuel  consumption  [2-4],  higher  reaction  zone  vol¬ 
ume  [5,6],  lower  combustion  noise  [7],  more  uniform  temperature 
field,  smaller  oscillations  of  temperature  in  the  reaction  zone  [7,8], 
sometimes  colorless  oxidation  of  fuel  [7,9],  lower  heat  release  rate 
and  lower  Damkohler  number  in  comparison  with  the  ordinary 
combustion  regime  [8,10,11],  and  the  influential  role  of  molecular 
diffusion  in  the  MILD  combustion  regime  [12]  . 

Dally  et  al.  [13]  reported  CO  measurements  in  an  experimental 
burner  (JHC)  which  emulates  the  MILD  combustion  mode.  CO  pro¬ 
files  at  different  02  concentration  show  a  reduction  at  lower  02  lev¬ 
els.  They  mentioned  that  there  is  a  different  mechanism  and 
chemical  pathway  in  MILD  flames  in  comparison  with  the  standard 
atmospheric  flames.  Numerical  modeling  of  the  JHC  burner  using  a 
RANS  based  code  has  been  conducted  by  Christo  and  Dally  [14],  in 
2005.  They  predicted  the  trends  of  measurements,  successfully. 
Christo  reported  that  the  CO  profile  has  an  unexpected  behavior 
by  02  concentration  reduction.  Because  higher  02  concentration 
in  the  high  preheated  co-flow  stream  is  expected  to  yield  higher 
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Fig.  1.  Computational  domain. 
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Fig.  2.  (a)  Methane  oxidation  pathways  for  oxidizer  of  23%  02  +  77%  N2  (by  mass)  in  a  WSR  analysis  (Case  1  in  Table  2)  at  T  =  1345  K  and  P  =  1  atm  for  a  0.1 -s  residence  time 
and  stoichiometric  condition.  Reaction  numbers  refer  to  Table  1,  while  reaction  rates  are  shown  in  parentheses  in  gmol/cm3  s.  (b  and  c)  Comparison  of  predicted  mole  fraction 
of  C02,  H20,  H,  OH,  O,  CH4,  CH3i  CH2  and  CO  for  methane  oxidation  by  oxidizer  of  23%  02  +  77%  N2  (by  mass)  in  a  WSR  (Cases  2,  3  and  4  in  Table  2)  by  four  chemical 
mechanisms  of  DRM-22,  GRI2.11,  C3_41,  and  NC5_44  under  conditions  of  P  =  1  atm  and  a  0.1 -s  residence  time  and  stoichiometric  ratio. 
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Table  1 

DRM-22  chemical  mechanism. 


No. 

Reaction 

No. 

Reaction 

No. 

Reaction 

1 

O  +  H  +  M^OH  +  M 

36 

H  +  HC0(+M)  o  CH20(+M) 

71 

H02  +  CO  ^  OH  +  C02 

2 

o  +  h2^h  +  oh 

37 

H  +  HCO  o  H2  +  CO 

72 

H02  +  CH20  ^  HCO  +  H202 

3 

0  +  H02  OH  +  02 

38 

H  +  CH20(+M)  o  CH30(+M) 

73 

CH2  +  02  ^  OH  +  HCO 

4 

0  +  CH2  ^  H  +  HCO 

39 

H  +  CH20  o  HCO  +  H2 

74 

CH2  +  H2  ^  H  +  CH3 

5 

0  +  CH2(S)^H  +  HC0 

40 

H  +  CH30  o  OH  +  CH3 

75 

2CH2  ^  H2  +  C2H2 

6 

0  +  CH3  ^  H  +  CH20 

41 

H  +  C2H2(+M)oC2H3(+M) 

76 

ch2  +  ch3  ^  h  +  c2h4 

7 

0  +  CH4  ^  OH  +  CH3 

42 

H  +  C2H3(+M)  o  C2H4(+M) 

77 

CH2  +  CH4  o  2CH3 

8 

0  +  CO  +  M  <3-  C02  +  M 

43 

H  +  C2H3  o  H2  +  C2H2 

78 

CH2(S)  +  n2  «•  ch2  +  n2 

9 

0  +  HCO  ^  OH  +  CO 

44 

H  +  C2H4(+M)  o  C2H5(+M) 

79 

CH2(S)  +  AR  44>  CH2  +  AR 

10 

0  +  HCO  ^  H  +  C02 

45 

H  +  C2H4  o  C2H3  +  H2 

80 

CH2(S)  +  02^H  +  0H  +  C0 

11 

0  +  CH20  ^  OH  +  HCO 

46 

H  +  C2H5(+M)  o  C2H6(+M) 

81 

CH2(S)  +  02^C0  +  H20 

12 

0  +  C2H2^CH2(S)  +  C0 

47 

H  +  C2H6  o  C2H5  +  H2 

82 

CH2(S)  +  H2^CH3  +  H 

13 

0  +  C2H2  o  CO  +  ch2 

48 

H2  +  C0(+M)oCH20(+M) 

83 

CH2(S)  +  H20^CH2  +  H20 

14 

0  +  C2H4  o  CH3  +  HCO 

49 

OH  +  H2  o  H  +  H20 

84 

CH2(S)  +  CH3^H  +  C2H4 

15 

0  +  c2h5  ^  ch3  +  ch2o 

50 

20H(+M)  o  H202(+M) 

85 

CH2(S)  +  CH4^2CH3 

16 

0  +  c2h6  ^  oh  +  c2h5 

51 

20H  o  0  +  H20 

86 

CH2(S)  +  CO^CH2  +  CO 

17 

02  +  CO  -<=>■  0  +  co2 

52 

OH  +  H02  o  02  +  H20 

87 

CH2(S)  +  C02^CH2  +  C02 

18 

02  +  CH20  o  H02  +  HCO 

53 

OH  +  H202  o  H02  +  H20 

88 

CH2(S)  +  C02^C0  +  CH20 

19 

H  +  02  +  M  o  H02  +  M 

54 

OH  +  CH2  o  H  +  CH20 

89 

ch3  +  02  ^  0  +  ch3o 

20 

H  +  202  <o  H02  +  02 

55 

0H  +  CH2(S)oH  +  CH20 

90 

ch3  +  o2^oh  +  ch2o 

21 

H  +  02  +  H20  o  H02  +  H20 

56 

OH  +  CH3  o  CH2  +  H20 

91 

ch3  +  H202  o  ho2  +  ch4 

22 

H  +  02  +  N2  <o  H02  +  N2 

57 

0H  +  CH3oCH2(S)  +  H20 

92 

2CH3(+M)^C2H6(+M) 

23 

H  +  02  +  AR  <o  H02  +  AR 

58 

OH  +  CH4  o  CH3  +  H20 

93 

2CH3  o  H  +  C2H5 

24 

H  +  02  o  0  +  OH 

59 

OH  +  CO  o  H  +  C02 

94 

CH3  +  HCO  ^  CH4  +  CO 

25 

2H  +  M  o  H2  +  M 

60 

OH  +  HCO  o  H20  +  CO 

95 
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26 

2H  +  H2  o  2H2 

61 

OH  +  CH20  o  HCO  +  H20 

96 

ch3  +  c2h4  ^  c2h3  +  ch4 

27 

2H  +  H20  o  H2  +  H20 

62 

OH  +  C2H2  o  CH3  +  CO 

97 

ch3  +  C2H6  o  c2h5  +  ch4 

28 

2H  +  C02  o  H2  +  C02 

63 

OH  +  C2H3  ^  H20  +  C2H2 

98 

HCO  +  H20  ^  H  +  CO  +  H20 

29 

H  +  OH  +  M  o  H20  +  M 

64 
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99 
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Fig.  3.  Comparison  of  radial  profiles  of  numerical  prediction  and  experimental  measurements  for  two  oxidizer  concentrations  of  3%  and  9%  02  at  Z  =  30  mm  (Cases  5  and  6  in 
Table  2). 


conversion  of  CO  to  C02.  They  reported  such  an  observation  under 
laminar  MILD  flame  and  guessed  that  the  pathway  of  CO  conver¬ 


sion  to  C02  could  change  under  MILD  condition,  de  Joannon  et  al. 
[15],  using  a  well-stirred  reactor  (WSR)  configuration,  reported 
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Table  2 

Numerical  experiments. 


No. 

Fuel 

Hot  oxidizer  composition 

CFD 

WSR 

Chemical  mechanism 

Conditions  of  WSR  (P^/I^.E4) 

1 

100%  ch4 

23%  02  +  77%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

2 

100%  ch4 

23%  02  +  77%  N2 

* 

GRI2.11 

(1,0.1,1345,1) 

3 

100%  ch4 

23%  02  +  77%  N2 

* 

C3-41 

(1,0.1,1345,1) 

4 

100%  ch4 

23%  02  +  77%  N2 

* 

NC5-44 

(1,0.1,1345,1) 

5 

20%  H2  +  80%  CH4 

9%  02  +  6.5%  H20  +  5.5%  C02  +  79%  N2 

* 

DRM-22 

- 

6 

20%  H2  +  80%  CH4 

3%  02  +  6.5%  H20  +  5.5%  C02  +  85%  N2 

* 

DRM-22 

- 

7 

100%  ch4 

9%  02  +  91%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

8 

100%  ch4 

3%  02  +  97%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

9 

20%  H2  +  80%  CH4 

9%  02  +  6.5%  H20  +  5.5%  C02  +  79%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

10 

20%  H2  +  80%  CH4 

3%  02  +  6.5%  H20  +  5.5%  C02  +  85%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

11 

20%  H2  +  80%  CH4 

9%  02  +  91%  N2 

* 

DRM-22 

(1.0.1.1345.0.8-1.2) 

12 

100%  ch4 

[3-23]%  02+  [97-77]%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

13 

20%  H2  +  80%  CH4 

[3-23]%  02+  [97-77]%  N2 

* 

DRM-22 

(1,0.1,1345,1) 

14 

100%  ch4 

100%  02 

* 

DRM-22 

(1.0.1.1200-2200.1) 

15 

20%  H2  +  80%  CH4 

100%  02 

* 

DRM-22 

(1.0.1.1200-2200.1) 

16 

100%  ch4 

3%  02  +  97% 

* 

DRM-22 

(1.0.1.1200-2200.1) 

17 

20%  H2  +  80%  CH4 

3%  02  +  97% 

* 

DRM-22 

(1.0.1.1200-2200.1) 

18 

20%  H2  +  80%  CH4 

23%  02  +  77%  N2 

* 

DRM-22 

(1.0.01-0.28.1345.1) 

19 

20%  H2  +  80%  CH4 

9%  02  +  91%  N2 

* 

DRM-22 

(1.0.01-0.28.1345.1) 

20 

20%  H2  +  80%  CH4 

3%  02  +  97%  N2 

* 

DRM-22 

(1.0.01-0.28.1345.1) 

1:  Pressure  in  atm;  2:  Residence  time  in  sec;  3:  Temperature  in  Kelvin;  4:  Equivalence  Ratio. 


that  the  oxidation  and  pyrolitic  regimes  are  the  working  conditions 
for  diluted  combustion.  The  former  is  suitable  for  reburning  tech¬ 
nology  which  results  in  CO  and  H20  as  main  reaction  products. 
The  latter  could  be  considered  for  a  reasonable  explanation  of 
flameless  combustion,  which  leads  to  large  production  of  CO  and 
H2.  Extremely  reduction  of  NOx,  CO,  and  UHC  emissions  under 
flameless  condition  is  also  reported  by  Li  et  al.  [16],  in  a  multiple 
jets  premixed  gas  turbine  combustor  and  a  swirling  jet  combustor. 
Recirculation  in  MILD  combustion  includes  C02  and  steam  in  reac¬ 
tion  zone.  Glarborg  and  Bentzen  [17]  and  Sabia  et  al.  [18]  showed 
that  carbon  dioxide  and  steam  can  modify  the  chemical  pathways 
involved  in  oxidation  processes.  Chinnici  et  al.  [19]  reported  the 
presence  of  temperature  oscillations  in  the  MILD  combustion  sys¬ 
tem  due  to  carbon  dioxide  in  the  reaction  zone,  de  Joannon  et  al. 
[20],  in  a  Jet  Stirred  Flow  Reactor,  investigated  the  temperature 
oscillation  in  methane  MILD  combustion  experimentally  and 
numerically.  They  reported  that  the  competition  between  oxida¬ 
tion  and  recombination  channels  is  stressed  under  the  high  dilu¬ 
tion  and  preheated  condition.  This  competition  is  responsible  for 
temperature  oscillation  modulation  in  stoichiometric  and  lean 
conditions.  Dally  et  al.  [21]  extended  the  applicability  of  MILD 
combustion  to  solid  biomass  fuel  combustion  in  a  laboratory  bur¬ 
ner.  Their  measurement  of  combustion  chamber  exhaust  gas  emis¬ 
sions  showed  the  reduction  of  NOx  emission  and  increment  of  CO 
by  transition  from  conventional  to  MILD  combustion  using  natural 
gas.  A  MILD  combustion  burner  fed  with  H2  and  H2/C02/H20  mix¬ 
tures  studied  by  Shabanian  et  al.  [22]  both  numerically  and  exper¬ 
imentally.  They  reported  the  values  of  NOx  and  CO  within  the 
exhaust. 

Methane  oxidation  to  CO  and  C02  under  high  and  low  temper¬ 
ature  has  been  studied  in  different  Refs.  [23,24],  but  CO  and  C02 
formation  under  MILD  condition  could  be  still  studied  in  more  de¬ 
tails.  Indeed  this  study  is  in  continuation  of  our  previous  reports  on 
MILD  combustion  [12,25,26],  and  its  main  objective  is  to  illustrate 
the  effect  of  oxygen  concentration  on  methane  oxidation  to  CO  and 
C02  under  MILD  combustion  regime.  Although  the  JHC  burner  con¬ 
figuration  has  shown  its  strength  to  emulate  MILD  combustion  un¬ 
der  controlled  condition  and  has  already  been  used  in  some  studies 
[27-30]  to  investigate  the  MILD  combustion  characteristics,  it  has 
been  less  mentioned  regarding  the  main  objective  of  present  work. 
In  this  framework,  a  CFD  analysis  of  JHC  burner  using  the  RANS 


equations  is  accompanied  by  zero  dimensional  analysis  of  a  WSR 
to  explore  the  carbon  monoxide  and  carbon  dioxide  formation  un¬ 
der  MILD  condition  for  different  oxidizer  02  levels  and  two  fuels  of 
methane  and  methane-hydrogen  mixture. 

2.  Numerical  method 

In  CFD  part,  the  JHC  burner  is  used  for  modeling.  The  numerical 
method  is  the  same  as  the  method  used  in  the  authors’  previous  re¬ 
ports  [12,25,26].  The  JHC  burner  is  a  co-flow  burner  which  is 
mounted  in  a  wind  tunnel  in  which  fuel  jet  issues  into  a  co-flow 
preheated  and  diluted  air.  Fuel  is  mixture  of  20%  H2  and  80%  CH4 
in  mass  basis.  The  fuel  jet  Reynolds  number  is  around  10,000  and 
the  velocity  of  hot  co-flow  air  is  3.2  m/s.  The  temperature  of  fuel, 
hot  co-flow  and  wind  tunnel  air  are  305,  1300,  and  305  K,  respec¬ 
tively.  Modeling  is  performed  using  the  2D-axisyymetric  RANS 
equations  and  modified  k-s  equations.  The  cei  coefficient  in  eddy 
dissipation  equation  is  changed  from  1.44  to  1.6,  according  to  the 
results  of  Christo  and  Dally  [14],  in  order  to  improve  the  standard 
k-s  model.  Velocity  profiles  at  the  inlets  are  estimated  from  non¬ 
reacting  flow  field  modeling  inside  the  burner.  According  to  the  re¬ 
sults  of  previous  reports  [12,14,25,31],  modeling  of  JHC  burner  is 
less  sensitive  to  the  turbulence  intensity  at  hot  co-flow  and  wind 
tunnel  inlets,  although  the  turbulence  intensity  at  the  fuel  inlet 
is  important  and  is  set  to  7%.  Patankar  SIMPLER  algorithm,  with  a 
third  order  Quick  discretization  method,  are  used  to  solve  the 
equations  [32].  Using  the  first  Fields  law  and  the  binary  diffusion 
coefficients,  found  from  molecular  kinetic  theory  [33],  differential 
diffusion  method  is  considered  but  the  radiation  effect  is  ne¬ 
glected.  Insignificant  effect  of  radiation  on  JHC  flame  modeling 
has  been  reported  by  Christo  and  Dally  in  2005  [14].  The  turbu¬ 
lence  and  chemistry  interaction  is  modeled  by  EDC  model.  The 
EDC  model  [34],  which  is  an  extension  of  the  eddy-dissipation 
model  (ED)  to  include  detailed  chemical  mechanisms  in  turbulent 
flows,  uses  finite  rate  model  in  conjunction  with  the  flow  held  tur¬ 
bulence  characteristics.  The  model  assumes  that  reaction  take 
places  in  the  small  turbulent  structures,  namely  the  fine  scales. 
The  volume  fraction  of  the  fine  scales  is  modeled  using  how  held 
turbulence  characteristics,  in  which  each  structure  is  a  Constant- 
Pressure-Reactor  (CPR)  and  homogeneous.  Species  are  assumed 
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Fig.  4.  Methane  oxidation  pathways  for  oxidizers  of  9%  02  +  91%  N2  (by  mass)  in  a  WSR  analysis  (Cases  7  in  Table  2)  at  T=  1345  K  and  P  =  1  atm  for  a  0.1 -s  residence  time  and 
stoichiometric  condition.  Reaction  numbers  refer  to  Table  1,  while  reaction  rates  are  shown  in  parentheses  in  gmol/cm3  s. 


to  react  in  the  fine  structures  over  a  time  scale  according  to  Arrhe¬ 
nius  finite  rate  law  and  detailed  chemical  mechanisms.  This  model 
is  used  extensively  in  MILD  combustion  modeling  [14,30,35]. 
Moreover,  its  performance  under  this  new  combustion  regime 
has  been  studied  in  detail  by  De  et  al.  [36].  Daily’s  burner  equiva¬ 
lent  CFD  computational  domain  is  shown  in  Fig.  1. 

Characteristics  of  the  MILD  combustion  regime  such  as  low 
reaction  rates  and  comparable  reaction  and  turbulence  time  scales 
[8,11]  seem  to  let  us  suppose  the  reaction  zone  composition  as  a 
outlet  of  Well-stirred  Reactor  (WSR)  [37].  The  characteristic  of 
the  WSR  entails  the  high  level  mixing  (i.e.  the  more  homogeneous 
mixture)  and  also  a  long  average  residence  time  inside  the  reactor 
(i.e.  finite  rate  reactions).  Linking  between  MILD  combustion  and  a 
WSR  characteristics  has  been  separately  reported  by  Plessing  et  al. 
[10]  and  also  Weber  et  al.  [38].  Attempt  of  De  Joannon  et  al.  [15]  for 
schematization  of  flameless  combustion  regime  has  been  made 
using  a  WSR  configuration.  They  concluded  that  the  fairly  similar¬ 
ity  between  concepts  of  MILD  and  WSR  are  more  related  to  oxida¬ 
tion  in  reburning  technology  in  terms  of  temperature  and 
composition.  On  the  other  hand  fine  structures,  in  the  EDC  model, 
are  assumed  as  isobaric,  adiabatic  and  perfectly  stirred  reactors. 
From  this  point  of  view,  it  is  possible  to  find  some  similarities  be¬ 
tween  characteristics  of  the  EDC  model  and  the  WSR  qualitatively. 
Therefore  the  WSR  analysis  could  be  a  useful  opportunity  to  be  a 
preliminary  analysis  to  catch  a  glimpse  of  methane  oxidation  pro¬ 
cess  under  low  oxygen  and  preheated  conditions  before  studying 
the  main  features  of  MILD  regime  by  reacting  flow  analysis,  here. 
The  WSR  analysis  is  performed  by  solving  the  species  conserva¬ 
tions  equations  under  constant  reactor  temperature  over  a  resi¬ 
dence  time  using  an  in-house  Fortran  code  which  includes  some 
subroutines  of  the  CFIEMKIN  open  source  code  [39]. 


3.  Chemical  mechanism 

To  study  CO  and  C02  formation,  following  the  reaction  path¬ 
ways  of  methane  oxidation  could  be  a  helpful  and  useful  approach. 
For  such  an  approach,  choosing  an  appropriate  and  accurate  chem¬ 
ical  mechanism  is  necessary  at  the  first  of  present  work.  To  this 


end,  there  have  been  many  detailed  chemical  mechanisms  devel¬ 
oped  to  describe  methane  oxidation.  They  consist  of  a  large  num¬ 
ber  of  species  and  reactions  which  using  them  in  reacting  flow 
modeling  entails  a  high  computational  cost.  Due  to  considerations 
of  our  computational  facilities,  choosing  a  costly  appropriate 
chemical  mechanism,  which  to  be  reasonably  able  to  describe 
the  main  features  of  methane  oxidation  under  MILD  condition,  is 
absolutely  vital  in  the  present  work.  Previous  reports  of  the  author 
and  the  other  different  reports  [25,26,40],  showed  an  acceptable 
accuracy  of  reduced  chemical  mechanism  of  DRM-22  [41  ]  in  com¬ 
parison  with  the  full  chemical  mechanism  of  GRI2.11. 

The  DRM-19  and  DRM-22  [41]  are  reduced  versions  of  the 
GRI1.2  [42],  which  their  performance  in  capturing  some  features 
of  MILD  combustion  has  been  identified  before.  Parente  et  al. 
[40]  reported  that  there  is  a  good  agreement  between  the  results 
of  the  DRM-19  and  the  GRI3.0  [43]  for  the  modeling  of  MILD  com¬ 
bustion.  On  the  other  hand,  under  MILD  condition,  in  addition  to 
the  uniformity  of  temperature  field,  extinction  and  re-ignition  oc¬ 
cur  broadly  in  the  reaction  zone  as  a  result  of  low  oxygen  concen¬ 
tration  and  high  temperature  preheating.  Therefore,  a  selected 
chemical  mechanism  should  be  able  to  capture  such  as  phenom¬ 
ena.  According  to  the  report  by  Kazakov  and  Frenklach  [41],  the 
DRM-22  shows  a  better  performance  than  DRM-19  in  predicting 
ignition  delay  time  and  laminar  flame  speed  in  atmospheric  pres¬ 
sure.  Therefore  it  seems  that  DRM-22  could  be  a  good  candidate 
for  representing  the  chemical  mechanism  in  the  MILD  combustion 
modeling.  In  this  way,  considering  some  more  detailed  analysis 
could  be  beneficial  to  increase  confidence  in  the  reasonable  accu¬ 
racy  of  DRM-22  in  the  prediction  of  methane  oxidation  pathways 
especially  relating  to  minor  species,  for  the  purposes  of  present 
work  here.  The  latter  concern  has  been  investigated  by  some  more 
detailed  analysis  as  follows. 

Methane,  which  exhibits  some  unique  combustion  characteris¬ 
tics  like  high  ignition  temperature  and  low  flame  speed,  is  widely 
researched  into  the  comprehensive  chemical  mechanism  [23]. 
Moreover  much  effort  has  been  directed  to  attain  simplified  mech¬ 
anisms  from  larger  comprehensive  mechanisms  for  different  com¬ 
bustion  modes.  At  first,  as  the  main  aim  of  the  present  work,  an 
attempt  has  been  made  to  choose  a  suitable  chemical  mechanism 
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Fig.  5.  Oxidation  pathways  of  80%  CH4  +  20%  H2  (by  mass)  for  two  oxidizers  of  (a)  9%  02  +  5.5%  C02  +  6.5%  H20  +  79%  N2,  (b)  3%  02  +  5.5%  C02  +  6.5%  H20  +  85%  N2  (by  mass)  in 
a  WSR  analysis  (Cases  9  and  10  in  Table  2)  at  T=  1345  K  and  P=  1  atm  for  a  0.1-s  residence  time  and  stoichiometric  condition.  Reaction  numbers  refer  to  Table  1,  while 
reaction  rates  are  shown  in  parentheses  in  gmol/cm3  s. 


with  regard  to  characteristic  of  MILD  combustion  and  limitation  of 
our  computational  facilities.  In  the  MILD  region  of  JHC  burner,  tem¬ 
perature  field  has  almost  a  low  gradient  around  1500  K  and  in  a 
whole  range  of  Low-Temperature  combustion  analysis  of  methane 
with  air  (i.e.  <1500  K)  which  is  explored  by  Glarborg  et  al.  [44]. 
Glarborg  and  coworkers  [23,44]  illustrated  the  major  chemical 
pathways  in  the  conversion  of  methane  to  C02  in  a  well-stirred 


reactor  at  temperature  of  1345  K  and  pressure  of  1  atm  for  a 
0.1  s  residence  time  using  the  GRI2.11  chemical  mechanism.  Under 
explored  condition,  they  reported  some  important  features  of  Low- 
Temperature  reaction  pathways  of  methane  like  a  strong  recombi¬ 
nation  of  CH3  back  to  CH4,  conversion  of  CH3  to  CH20  through  the 
intermediate  production  of  CH30,  and  producing  a  higher  hydro¬ 
carbon  than  the  original  reactant  methane  (i.e.  Ethane  and 
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Fig.  6.  Variation  of  CH3,  CO,  and  C2H6  rates  of  production  by  mentioned  reactions  in 
figure  due  to  changing  of  equivalence  ratio  for  WSR  calculation  at  T=1345K, 
P=  1  atm  and  a  residence  time  of  0.1  s  (Case  11  in  Table  2). 


Acetylene).  Similar  to  the  analysis  of  Glarborg  et  al.  [23,44],  using 
the  zero  dimensional  analysis  and  DRM-22  chemical  mechanism, 
pathways  of  pure  methane  (100%  CH4)  oxidation  to  CO  and  C02 
are  shown  in  Fig.  2a.  It  is  for  oxidizer  of  23%  02  +  77%  N2  (by  mass), 
constant  temperature  of  1345  K,  residence  time  of  0.1  s,  atmo¬ 
spheric  pressure,  and  stoichiometric  condition.  In  Fig.  2a,  each  ar¬ 
row  indicates  a  conversion  of  species  at  its  tail  to  the  species  at  the 
head.  The  mentioned  conversion  is  addressed  by  elementary  reac¬ 
tions  which  are  indicated  along  the  length  of  the  arrow  by  the 
numbers  of  reactions  in  the  DRM-22  mechanism  (Table  1)  and 
the  additional  reactant  species  in  the  mentioned  reaction,  while 
the  rate  of  destruction  of  reactant  is  quantified  by  the  parenthetical 
numerical  value.  Furthermore  the  width  and  color  of  arrows  gives 
a  visual  indication  of  the  relative  importance  of  each  pathway.  The 
elementary  reactions  with  rates  larger  than  1.0e-9  (gmol/cm3  s) 
are  considered  for  pathway  extractions  and  the  other  reactions 
are  filtered.  It  can  be  concluded  from  Fig.  2a  that  the  main  features 
of  methane  combustion  in  the  condition  of  Low-Temperature, 
which  are  mentioned  above,  have  been  captured  by  DRM-22  with 
a  lower  computational  cost  and  almost  same  accuracy  as  the  more 
complicated  and  comprehensive  GRI2.11  mechanism. 

Moreover  Curran  et  al.  [45]  has  mentioned  the  necessity  of 
including  the  CH302  radical  species  and  reactions  in  a  reaction 
mechanism  in  order  to  correctly  simulate  the  methane  oxidation 
chemistry  under  high-pressure  and  intermediate  temperature  con¬ 
ditions.  The  CH302  radical  species  and  reactions  are  not  included  in 
GRI-MECF1  family.  Flowever,  the  conditions  of  present  work  are  not 
in  the  range  of  high  pressure  but  it  is  in  the  range  of  intermediate 
temperature  as  mentioned  above.  To  consider  the  latter  point  con¬ 
cerning  the  performance  of  DRM-22  chemical  mechanism  over  the 
condition  of  intermediate  temperature,  some  more  studies  are 
planned  as  follows.  The  Low-Temperature  behavior  of  DRM-22  is 
studied  in  more  detail  by  comparing  its  results  with  three  more 
comprehensive  detailed  chemical  mechanisms  in  Fig.  2b  and  c.  De¬ 
picted  in  Fig.  2b  and  c  are  CO,  C02,  H,  OH,  O,  CH4,  CH3,  CH2  and  H20 
predicted  profiles  result  of  the  WSR  analysis  for  fuel  of  methane 
and  oxidizer  of  23%  02  +  77%  N2  (by  mass)  in  the  temperature 
range  of  1300-2000  K,  residence  time  of  0.1  s,  atmosphere  pres¬ 
sure,  and  stoichiometric  condition. 

These  detailed  chemical  mechanisms  are  GRI2.11,  C3_41,  and 
NC5_44.  C3_41  and  NC5_44  are  detailed  chemical  mechanism 
which  are  introduced  by  Curran  et  al.  [46]  for  natural  gas.  C3_41 


is  a  complete  kinetic  mechanism  for  natural  gas  which  consists 
of  118  different  chemical  species  and  663  elementary  reactions. 
It  includes  three  sub-mechanisms  for  hydrogen,  methane/ethane, 
and  C3.  The  hydrogen  sub-mechanism  is  based  on  what  was  vali¬ 
dated  by  O’Conaire  et  al.  [47]  in  the  temperature  range  of  298- 
2700  K,  at  pressures  from  0.05  to  87  atm,  and  equivalence  ratios 
from  0.2  to  6.0.  The  kinetic  mechanism  employed  for  the  meth¬ 
ane/ethane  system  is  based  on  that  published  by  Fischer  et  al. 
[48].  The  C3  sub-mechanism  is  based  on  the  work  of  Curran 
et  al.  [49,50]. 

NC5_44  is  a  detailed  chemical  mechanism  for  natural  gas 
including  the  higher-order  hydrocarbon  (for  example  ethane 
(C2H6),  propane  (C3H8),  butane  (C4Hi0),  etc.).  It  is  developed  by 
Curran  and  co-workers  [51-53].  This  mechanism  which  consists 
289  species  and  1580  reactions  has  produced  quite  well  results 
for  different  mixtures  at  a  wide  range  of  pressure  intervals 
(1-50  atm)  over  the  temperature  range  740-1550  K.  It  can  be 
understood  from  Fig.  2b  and  c  that  the  results  of  four  chemical 
mechanisms  are  in  a  reasonable  agreement  with  each  other  in 
the  temperature  range  of  1300-2000  K,  although  some  species  like 
CH2  and  C2H6  are  underpredicted  using  the  GRI-MECH  family  in 
comparison  with  more  complicated  mechanisms  of  NC5_44. 
C3_41  and  NC5_44. 

In  general,  the  DRM  mechanism  includes  almost  all  key  path¬ 
ways  through  which  methane  oxidation  is  thought  to  occur.  There¬ 
fore  due  to  lower  computational  cost,  DRM-22  is  used  as  the  main 
chemical  mechanism  in  the  present  work. 


4.  Results 

4.1.  Validation 

A  structured  grid  with  about  39,000  cells  is  selected  for  calcula¬ 
tions  after  a  grid  study.  CFD  results  have  been  verified  by  compar¬ 
ing  with  measurements  of  Dally  et  al.  [13]  in  Fig.  3.  Good 
agreement  between  numerical  and  experimental  data  for  temper¬ 
ature  and  main  species  of  H20  and  C02  is  a  good  sign  of  reasonable 
accuracy  of  modeling  of  flame  structure.  Furthermore  acceptable 
prediction  of  minor  species  like  OH  and  CO  could  be  an  evidence 
of  proper  performance  of  intermediate  reactions  of  chemical  mech¬ 
anism.  Twenty  numerical  test  cases  (Table  2)  are  defined  to  con¬ 
duct  the  research  toward  the  main  objective  of  the  present  work. 
They  include  two  test  cases  for  CFD  analysis  to  investigate  the  ef¬ 
fect  of  oxygen  concentration  on  reaction  zone  and  18  test  cases 
for  WSR  analysis  to  explore  the  effect  of  fuel  hydrogen  content, 
oxidizer  oxygen  concentration,  mixture  inlet  temperature,  and 
equivalence  ratio  on  methane  oxidation  process  to  CO  and  C02  un¬ 
der  MILD  condition. 


4.2.  The  well  stirred  reactor  (WSR )  analysis 

The  reasonable  result  of  Low-Temperature  analysis  of  methane 
oxidation,  using  DRM-22,  was  discussed  later  in  Fig.  2.  Importance 
of  ethane  in  methane  oxidation  pathways  under  the  calculated 
condition  is  the  most  important  finding  of  this  figure.  Ethane  route 
converts  the  methyl  to  CO  through  the  acetylene  while  the  main 
methyl  oxidation  route  is  CH3  -►  CH20  -►  HCO  -►  CO  -►  C02.  These 
pathways  were  also  reported  by  Glarborg  et  al.  [44].  In  continua¬ 
tion,  the  well-stirred  reactor  analysis  is  considered  to  study  the 
combustion  under  MILD  condition.  Such  an  attempt  provides  a 
general  overview  of  reaction  pathways  of  methane  oxidation  under 
MILD  condition.  At  first,  the  condition  of  well-stirred  reactor  calcu¬ 
lation  must  be  determined.  This  condition  could  be  defined  using 
the  results  of  CFD  modeling  of  JHC  burner.  Although  the  reacting 
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Fig.  7.  (a)  First  sixth  important  CO  containing  reactions  which  are  selected  by  sorting  CO  containing  reactions  according  to  the  absolute  rate  of  reactions  and  (b)  First  sixth 
important  C02  containing  reactions  which  are  selected  by  sorting  C02  containing  reactions  according  to  the  absolute  rate  of  reactions,  in  CFD  modeling  at  each  radial  distance 
for  Z  =  30  mm  and  9%  02  configuration  (Case  5  in  Table  2). 
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Fig.  8.  Reactive  flow  analysis  for  80%  CH4  +  20%  H2  (by  mass)  and  oxidizers  of  9%  02  +  5.5%  C02  +  6.5%  H20  +  79%  N2  (a)  MILD  region  in  OH  and  temperature  contours  and  (b) 
averaged  reaction  pathways  in  the  zone  of  part  a.  Reaction  numbers  refer  to  Table  1,  while  reaction  rates  are  shown  in  parentheses  in  gmol/cm3  s  (Case  5  in  Table  2). 


flow  characteristics  like  as  residence  time,  temperature,  and  fuel/ 
oxidizer  equivalence  ratio,  with  the  exception  of  pressure,  change 
locally  through  the  MILD  combustion  zone  (this  zone  is  deter¬ 
mined  and  discussed  in  Section  4.3),  it  has  been  tried  to  consider 
some  appropriate  values  as  a  whole.  Flow  residence  time  using 
the  formulation  of  EDC  model,  t*  =  cT(v/e)1/2,  changes  approxi¬ 
mately  between  0.02  and  1.8  s  through  the  MILD  zone.  As  an  aver¬ 
age,  the  residence  time  could  be  around  0.9  s.  In  the  t*  equation,  v 
is  kinematic  viscosity,  s  is  turbulent  dissipation  rate,  and 
cT  =  0.4082  is  the  time  scale  constant.  Furthermore  average  tem¬ 
perature  in  MILD  zone  is  around  the  values  ^  1500  K,  which  is  in 
the  temperature  range  of  Low-Temperature  oxidation  of  methane. 
To  have  a  consistency  between  conditions  of  WSR  calculations  of 
all  cases  here  and  also  noticing  that  the  aim  of  WSR  analysis  in 
present  work  is  more  qualitative  rather  quantitative,  WSR  condi¬ 


tions  are  set  to  constant  temperature  of  1345  K,  residence  time 
of  0.1  s,  atmospheric  pressure,  and  stoichiometric  condition  same 
as  the  setting  of  Case  1  in  Table  2. 

MILD  combustion  occurs  under  much  lower  02  concentration 
than  23%  02  (by  mass).  Therefore,  comparing  to  Case  1,  02  concen¬ 
tration  is  decreased  from  23%  to  9%  in  Fig.  4.  In  comparison  with 
atmospheric  condition  (Fig.  2),  rates  of  reactions  are  reduced  but 
ethane  still  has  a  key  role  in  methane  oxidation.  While  the  con¬ 
necting  reactions  between  ethane  route  and  main  CF14  oxidation 
pathway  (i.e.  CH3  CH20  HCO  -»  CO  -►  C02)  and  also  CFI3O 
containing  reactions  are  weakened.  On  the  other  hand,  pathway 
of  CH2  CH20  weakens  although  CFI2  HCO  intensifies.  Another 
important  point  is  that  the  rates  of  reactions  R59  and  R12-13  have 
decreased  by  79%  and  99%,  respectively,  due  to  reduction  in  inlet 
02  concentration.  That  means  a  lower  relative  reduction  of  main 


A.  Mardani  et  al.  /  Combustion  and  Flame  160  (2013)  1636-1649 


1645 


02  MoleFraction  in  Inlet  Air  of  WSR 


Fig.  9.  Effect  of  oxidizer  02  concentration  on  formation  of  CO  and  C02  for  two  fuels 
of  100%  CH4  and  80%  CH4  +  20%  H2  and  oxidizer  of  02/N2  under  stoichiometric 
condition  in  a  WSR  at  T=  1345  K  and  P=  1  atm  for  a  0.1-s  residence  time  (Cases  12 
and  13  in  Table  2). 
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Fig.  11.  Effect  of  WSR  temperature  on  formation  of  CO  and  C02  for  two  fuels  of 
100%  CH4  and  80%  CH4  +  20%  H2  and  oxidizer  of  3%  02  +  97%  N2  (by  mass)  under 
stoichiometric  condition  in  a  WSR  at  P  =  1  atm  for  a  0.1-s  residence  time  (Cases  16 
and  17  in  Table  2). 
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Fig.  10.  Effect  of  WSR  temperature  on  formation  of  CO  and  C02  for  two  fuels  of 
100%  CH4  and  80%  CH4  +  20%  H2  and  oxidizer  of  100%  02  under  stoichiometric 
condition  in  a  WSR  at  P  =  1  atm  for  a  0.1-s  residence  time  (Cases  14  and  15  in 
Table  2). 


Fig.  12.  Effect  of  residence  time  in  WSR  on  formation  of  CO  and  C02  for  fuel  of  80% 
CH4  +  20%  H2  and  three  oxidizers  of  3%  02  +  97%  N2>  9%  02  +  91%  N2f  23%  02  +  77%  N2 
(by  mass)  under  stoichiometric  condition  in  a  WSR  at  P  =  1  atm  (Cases  18, 19  and  20 
in  Table  2). 


C02  producer  reaction  (i.e.  R59)  in  respect  to  CO  producer  reactions 
R12  and  R13. 

More  decreasing  of  02  concentration  to  3%  02  under  the  other 
constant  conditions,  reduces  the  reaction  rates  significantly  to  low¬ 
er  than  1.0e-9  (gmol/cm3  s),  therefore  the  reaction  pathways  for 
3%  02  are  ignored  to  present  here.  This  (i.e.  the  ignorable  reaction 
rates  for  3%  02)  is  a  good  reason  for  adding  H2  to  methane  in  the 
JHC  burner  fuel  feeding.  Because  as  which  mentioned  in  Ref.  [12], 
hydrogen  addition  to  methane  is  necessary  for  sustaining  the  flame 
at  the  JHC  setup  for  3%  02  and  the  rate  of  reactions  become  ignor¬ 
able  below  2%  H2  (by  mass  fraction).  This  would  be  discussed  more 
as  follows.  In  the  next  step,  methane  is  enriched  by  20%  H2  and 
methane  oxidation  pathways  to  CO  and  C02  are  illustrated  for  9% 


02  in  Fig.  5a.  It  can  be  seen  that  the  production  of  intermediate  spe¬ 
cies  of  formaldehyde  (CH20),  which  could  be  considered  as  an  igni¬ 
tion  marker,  has  increased  as  a  result  of  intensified  pathways  of 
CH2  CH20  and  CH3  CH20.  Formation  of  CH2  increases  due  to 
direct  conversion  of  CH3  through  reaction  R56  and  also  route  of 
CH3  CH2(s)  CH2.  Moreover  an  alternate  route  from  CH3  to 
CH20  appears  through  the  intermediate  species  of  CH30.  It  shows 
that  ignition  process  is  more  powerful  for  blended  methane  with 
H2  in  comparison  with  pure  methane.  It  is  in  consistency  with 
the  result  of  Sabia  et  al.  [54]  on  the  higher  the  increase  of  the  sys¬ 
tem  reactivity  with  the  higher  fuel  inlet  hydrogen  content.  Indeed, 
a  comparison  between  Figs.  4  and  5a  reveals  that  hydrogen  modi¬ 
fies  the  evolution  of  the  methane  oxidation  process  through  the 
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Fig.  13.  (a)  First  sixth  important  CO  containing  reactions  which  are  selected  by  sorting  CO  containing  reactions  according  to  the  absolute  rate  of  reactions,  (b)  First  sixth 
important  C02  containing  reactions  which  are  selected  by  sorting  C02  containing  reactions  according  to  the  absolute  rate  of  reactions,  in  CFD  modeling  at  each  radial  distance 
for  Z  =  30  mm  and  3%  02  configuration  (Case  6  in  Table  2). 
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Fig.  14.  Comparison  of  reaction  rates  of  R13  and  R59  along  a  line  inside  the  reaction 
zone  (point  of  maximum  OH  mass  fraction  at  each  distance  from  the  fuel  nozzle) 
inside  the  MILD  combustion  zone  (Z<100mm,  according  to  Ref.  [13])  for  two 
oxidizer  of  3%  and  9%  (Cases  5  and  6  in  Table  2). 
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Fig.  15.  Comparison  of  radial  profiles  of  0,  OH,  and  H  radicals  at  Z  =  30  mm  for  two 
oxidizer  of  3%  and  9%  02  (Cases  5  and  6  in  Table  2). 


growth  in  reaction  rates  of  the  main  methane  oxidation  pathways 
to  C02  (i.e.  CH3  CH20  HCO  -►  CO  -►  C02)  and  also  a  decline  in 
reaction  rates  of  alternative  methane  oxidation  route  to  CO 
through  the  higher  hydrocarbon  than  the  original  reactant  meth¬ 
ane  (i.e.  CH3  -  C2H6  -  C2H5  -  C2H4  -  C2H3  -  C2H2  -  CO). 

In  continuation,  reactions  pathways  for  fuel  of  80%  CH4  +  20%  H2 
and  oxidizer  of  3%  02  +  97%  N2  are  shown  in  Fig.  5b.  It  can  be 
understood  from  the  figure  that,  in  spite  of  pure  methane  and  3% 
02,  rate  of  reactions  are  considerable  as  a  result  of  presence  of  H2 
in  the  fuel  mixture.  Moreover,  by  decreasing  of  02  from  9%  to  3%, 
rate  of  reactions  decreases  and  CH30  and  ethane  oxidation  path¬ 
ways  could  be  ignored.  Furthermore,  reaction  rates  of  all  signifi¬ 
cant  pathways  are  in  the  same  order  of  magnitude.  That  means 
for  3%  02,  all  shown  pathways  (Fig.  5b)  have  an  approximately 
same  importance.  Weakening  of  pathways  ended  by  CO  and  reac¬ 
tion  R59  could  be  a  sign  of  lower  CO  and  C02  production  at  3%  02. 


On  the  whole,  the  WSR  analysis  reveals  that  the  Low-Tempera¬ 
ture  methane  oxidation  pathways  should  be  thought  over  under 
diluted  and  preheated  condition  (around  1500  K).  This  finding, 
which  will  be  noticed  in  the  continuation  of  this  work,  is  a  result 
of  WSR  calculation  under  stoichiometric  condition.  On  the  other 
hand,  numerical  modeling  of  JHC  burner  has  shown  the  occurrence 
of  maximum  temperature  or  OH  mass  fraction  in  the  vicinity  of 
stoichiometric  zone  and  not  exactly  there.  Therefore  it  would  be 
worthwhile  to  investigate  the  effect  of  equivalence  ratio  on  the 
CO  and  C02  formation  routes  to  make  sure  about  preserving  overall 
trends  in  mentioned  finding.  In  other  word,  it  is  important  to  know 
whether  the  higher  hydrocarbon  oxidation  route  is  still  important 
for  the  mixture  different  from  stoichiometric  ratio.  To  follow  such  a 
question,  variations  of  C2H6  production  from  CH3  (by  R92  in  the 
start  of  ethane  route)  and  CO  production  from  C2H2  (by  R12  and 
R13  at  the  end  of  ethane  route)  versus  inlet  mixture  equivalence 
ratio  are  depicted  in  Fig.  6  for  condition  of  9%  02.  Furthermore, 
for  more  detail,  comparison  between  conversions  of  CH3  to  C2H6 
(by  R92),  to  CH20  (by  R6),  and  to  CH2(s)  (by  R57)  are  also  shown 
in  Fig.  6  for  different  inlet  mixture  equivalence  ratio.  It  can  be 
understood  from  the  figure  that  in  the  rich  region  (equivalence  ra¬ 
tio  higher  than  1.0)  the  methyl  conversion  to  C2H6  increases  and 
CH20  and  CH2(s)  formations  decrease  and  also  CO  formation  by 
ethane  route  increases.  They  are  evidence  that  the  higher  hydro¬ 
carbon  route  grows  in  importance  for  rich  condition.  On  the  other 
hand,  under  lean  condition,  related  reactions  to  ethane  route  have 
approximately  same  rates  in  comparison  with  stoichiometric  con¬ 
dition.  That  means  that  it  is  possible  to  diversify  the  predicted 
methane  oxidation  pathways  (to  CO  and  C02)  under  stoichiometric 
condition  to  a  little  leaner  condition.  Indeed,  stoichiometric  WSR 
analysis  could  be  considered  as  an  indicator  whether  the  higher 
hydrocarbon  oxidation  route  is  still  important  for  the  mixture  dif¬ 
ferent  from  stoichiometric  ratio  (albeit  in  the  vicinity  of  stoichiom¬ 
etric  condition)  and  also  globally  could  be  used  to  explain  points 
regarding  CO  and  C02  formation  in  MILD  regime  reaction  zone 
even  for  non-stoichiometric  condition. 

4.3.  Reacting  flow  analysis 

In  this  section  changes  of  CO  and  C02  concentration,  due  to  var¬ 
iation  of  02  concentration,  are  studied  at  30  mm  from  the  nozzle  of 
JHC  burner.  The  distance  of  30  mm  is  considered  as  a  representa¬ 
tive  sample  of  MILD  region  because  it  is  inside  of  MILD  region  of 
reaction  zone  (Z<  100  mm  in  JHC).  In  addition,  Dally  et  al.  [13] 
have  reported  experimental  measurements  there.  In  this  frame¬ 
work,  reactions  of  CO  and  C02  are  focused. 

In  DRM-22,  CO  including  reactions  are  R8,  R10,  R9,  R12,  R13, 
R17,  R37,  R48,  R59,  R60,  R62,  R71,  R81,  R86,  R88,  R94,  R98,  R99, 
R100  and  C02  including  reactions  are  R8,  R10,  R17,  R28,  R59, 
R71,  R87,  and  R88.  Depicted  in  Fig.  7a  and  b  respectively  are  sorted 
CO  and  C02  producer  and  consumer  reactions  according  to  their 
absolute  rates  magnitude  at  each  radial  location  of  Z  =  30  mm  for 
9%  02.  Moreover,  the  types  of  reactions  (i.e.  CO  or  C02  producer 
reaction)  are  identified  in  attached  tables  as  well  as  their  sequence. 
They  show  that  in  the  reaction  zone  ( R  =  7.1-8.464  mm),  the  most 
dominant  C02  containing  reactions  are  R59,  R88,  and  R8  and  also 
the  paramount  important  reactions  including  CO  species  are  R59, 
R12,  R13,  R99,  and  R98.  The  temperature  profile  at  Z  =  30mm 
(Fig.  3)  shows  that  the  averaged  temperature  in  reaction  zone  is 
around  1500  I<  and  therefore  methane  reaction  pathways  could 
be  considered  broadly  similar  to  Fig.  5a  in  which  the  reactor  inlet 
oxygen  concentration  and  temperature  are  9%  and  1345  K,  respec¬ 
tively.  In  tables  of  Fig.  7,  symbolizing  and  dividing  the  reactions  to 
producer  and  consumer  of  CO  and  C02  are  conducted  by  consider¬ 
ing  methane  oxidation  pathways  in  Fig.  5a.  The  reaction 
arrangement  in  Fig.  7a  shows  that  the  reactions  of  C2H2  CO 
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(i.e.  R12-13)  are  more  important  to  CO  formation  in  comparison 
with  reactions  of  HCO  CO  (i.e.  R98-99)  and  CH2(s)  CO  (i.e. 
R88).  That  means  that  in  this  location  of  JHC  flame  and  under  con¬ 
sidered  condition,  ethane  route  attaches  great  importance  to  meth¬ 
ane  oxidation.  This  finding  is  more  studied  by  extracting  the 
chemical  pathways  in  reacting  flow  analysis  using  the  averaged 
reaction  rates  in  MILD  combustion  region.  MILD  zone  is  considered 
as  a  region  with  Z  <  100  mm  and  OH  mass  fraction  of  higher  than 
one  percent  of  maximum  OH  mass  fraction  at  each  axial  distance 
from  the  fuel  nozzle(Z).  Measurements  of  Dally  et  al.  [13]  on  JHC 
burner  have  illustrated  that  the  mixing  with  fresh  tunnel  air  affects 
the  flame  above  100  mm  from  the  nozzle  and  the  MILD  combus¬ 
tion  regime  occurs  at  the  region  below  this  height.  MILD  zone  is 
illustrated  in  Fig.  8a  by  dashed  close  line  for  jet  Reynolds  number 
of  10,000  and  02  mass  fraction  of  9%.  Reaction  pathways,  based 
on  averaged  rates  of  reactions  in  MILD  zone  (Fig.  8a),  are  depicted 
in  Fig.  8b.  It  can  be  understood  that  the  reacting  flow  analysis  also 
shows  the  unquestionable  importance  of  the  ethane  route  in  the 
methane  oxidation  to  CO  and  C02.  In  other  words  there  is  a  consis¬ 
tency  between  results  of  WSR  and  CFD  on  importance  of  high 
hydrocarbon  pathways  in  the  methane  oxidation  under  MILD 
condition. 

A  comparison  between  Fig.  7a  and  b  shows  that  the  most 
important  reaction  for  C02  formation  (i.e.  R59)  has  a  smaller  rate 
in  comparison  with  reactions  R12  and  R13.  That  means  that  the 
CO  producer  reactions  R12  and  R13  have  higher  rates  than  the 
CO  consumer  reaction  (or  C02  producer  reaction)  R59.  Moreover, 
the  reactions  R59  and  R88  are  significant  among  C02  containing 
reactions  which  R59  is  more  important  than  the  others. 

Figure  3  and  both  measurements  of  Dally  et  al.  [13]  and  calcu¬ 
lations  of  Christo  and  Dally  [14]  show  an  increment  in  CO  and  C02 
formation  by  increasing  of  02  concentration.  As  explained  in  the 
introduction,  the  behavior  of  CO  is  fairly  strange.  To  explain  this 
phenomenon,  the  behavior  of  CO  and  C02  concentration  variation 
due  to  the  changes  of  02  concentration  are  studied  in  a  WSR  first 
and  after  that  again  CFD  results  are  considered  to  complete  the  dis¬ 
cussion  about  CO  behavior. 

Figure  9  shows  the  correlations  between  CO  and  C02  concentra¬ 
tions  and  inlet  02  concentration  in  a  WSR  for  two  different  fuel 
compositions  under  stoichiometric  condition,  residence  time  of 
0.1  s,  and  temperature  of  1345  K.  The  figure  illustrates  the  incre¬ 
ment  in  CO  and  C02  concentrations  at  higher  02  levels.  Further¬ 
more,  the  rates  of  CO  consumer  reactions  also  show  an  increase 
for  higher  02  concentration  (not  shown  here).  On  the  other  hand 
the  WSR  temperature  is  another  important  parameter  which  its  ef¬ 
fect  on  CO  and  C02  concentration  is  studied  in  Figs.  10  and  11.  De¬ 
picted  in  Fig.  10  are  the  effects  of  combustion  chamber 
temperature  on  CO  and  C02  concentration  for  two  fuels  of  100% 
CH4  and  80%  CH4  +  20%  H2  and  oxidizer  of  100%  02.  It  can  be  seen 
that  the  effect  of  temperature  on  CO  and  C02  concentration  is  dif¬ 
ferent  at  temperature  higher  and  lower  than  approximately 
1500  K.  That  means,  under  condition  of  T  <  1500  K  a  decrease  in 
combustion  chamber  temperature  increases  the  CO  concentration 
and  decreases  C02  concentration,  while  for  T>  1500  K  it  is  quite 
the  reverse.  That  is  probably  due  to  dissociation  of  C02  at  temper¬ 
atures  higher  than  1500  K.  Figure  11  is  the  same  as  Fig.  10  but  only 
for  oxidizer  containing  3%  02.  It  shows  that  at  lower  02  concentra¬ 
tion,  combustion  chamber  temperature  has  an  approximately 
identical  effect  on  CO  and  C02  concentration  and  the  concentration 
of  both  CO  and  C02  increase  in  higher  combustion  temperature. 
Therefore,  the  latter  analysis  reveals  that  the  observed  behavior 
(i.e.  CO  variation)  in  JHC  burner,  due  to  increase  in  02  level,  could 
be  considered  approximately  equivalent  to  the  effect  of  flame  tem¬ 
perature  at  Fig.  11.  The  range  of  temperature  in  the  MILD  zone  of 
JHC  burner  passes  around  1500 1<  with  increasing  oxidizer 
concentration. 


To  get  a  better  understanding  of  CO  and  C02  variation  by  the 
changes  in  02  concentration,  the  effect  of  residence  time  on  con¬ 
centration  of  CO  and  C02  for  three  02  levels  of  3%,  9%,  and  21% 
are  depicted  in  Fig.  12.  It  shows  that  the  concentration  of  CO  and 
C02  respectively  decreases  and  increases  with  an  increase  of  resi¬ 
dence  time.  Moreover,  the  CO  concentration  is  more  sensitive  to 
residence  time  in  comparison  with  C02  for  all  02  concentrations, 
which  this  sensitivity  increases  for  higher  02  level.  That  means  that 
fluid  dynamics,  which  has  a  direct  impact  on  residence  time,  is 
effective  on  CO  and  C02  production.  This  correlation  probably 
grows  in  importance  for  higher  02  levels.  Therefore  the  residence 
time  is  responsible  for  some  portion  of  CO  concentration  increment 
at  higher  02  levels  in  the  JHC  burner.  It  should  be  considered  that 
MILD  combustion  and  WSR  have  similar  physical  characteristics. 
Furthermore  there  is  a  consistency  between  the  residence  time  ef¬ 
fects  at  higher  02  levels  and  the  effect  of  02  concentration  on  tur¬ 
bulence  relaminarization  which  is  discussed  extensively  in  the 
other  report  of  the  authors  [26].  Damping  turbulence  eddies  and 
relaminarizing  of  flow  increases  at  higher  02  concentrations  be¬ 
cause  the  rate  of  heat  release  increases  due  to  increase  in  02  level. 

In  continuation,  the  variation  of  CO  and  C02  concentration  due 
to  the  change  of  02  concentration  are  studied  using  the  CFD  results. 
In  a  similar  approach  as  done  for  Figs.  7  and  13  is  made  for  3%  02 
and  Z  =  30  mm.  This  figure  illustrates  that  the  C02  producer  and  CO 
consumer  reaction  R59  has  a  higher  rate  than  the  CO  producer 
reactions  of  the  ethane  route  (i.e.  R12  and  R13),  while  the  reverse 
can  be  understood  from  Fig.  7a  for  9%  02.  That  means  the  reactions 
R12  and  R13  have  higher  rates  than  R59  for  9%  02  and  they  show  a 
contrary  order  for  3%  02.  So,  it  could  be  understood  that  the  behav¬ 
iors  of  CO  and  C02  concentrations  at  different  02  levels  are  due  to 
correlation  between  reactions  R12-13  and  R59.  In  other  words, 
higher  hydrocarbon  oxidation  (i.e.  ethane)  is  responsible  for  CO 
and  C02  concentration  behavior  by  02  concentration  variation.  This 
idea  is  more  studied  by  comparing  reaction  rats  R12-13  and  R59 
along  a  line  inside  of  the  reaction  zone  of  the  jet  (i.e.  different  Z). 
Such  a  comparison  for  two  02  levels  of  3%  and  9%  are  depicted  in 
Fig.  14.  This  figure  shows  that  by  changing  02  concentration,  the 
ratio  between  reaction  rates  of  R13  and  R59  changes.  Furthermore, 
a  comparison  between  Figs.  7a  and  13a  reveals  that  another  CO 
producer  reaction  (i.e.  R88)  is  limited  at  lower  02  percent,  and  it 
might  be  due  to  lower  CH2(s)  production  from  CH3.  To  understand 
how  inlet  oxygen  concentration  affects  the  methane  oxidation 
pathways,  the  prominent  radicals  should  be  considered  more  pre¬ 
cisely.  Radical  O  has  an  important  role  in  R12-13,  which  its  con¬ 
centration  decreases  at  lower  02  levels.  Other  important  radicals 
like  OH  and  H  also  reduce  at  lower  02  concentrations  (see 
Fig.  15).  These  radicals  are  the  main  reason  of  reactions  suppres¬ 
sion  at  lower  oxidizer  02  levels.  Because  in  the  first  step,  oxidation 
process  starts  with  an  attack  on  the  CH4  molecule  by  OH,  O,  and  H 
radicals  to  produce  the  methyl  radical  and  in  continuation  these 
radicals  play  a  key  role  in  the  proceeding  of  both  main  and  higher 
hydrocarbon  oxidation  pathways.  In  the  main  route,  the  methyl 
radical  forms  formaldehyde  (CH20)  by  combining  with  an  oxygen 
atom  or  with  an  OH  and  H  radicals  through  the  intermediate 
CH2;  the  formaldehyde,  in  turn,  is  attacked  by  OH,  H,  and  O  radicals 
to  produce  the  formyl  radical(HCO).  In  the  ethane  route,  OH  attacks 
to  C2H6  and  C2H4  to  form  C2H5  and  C2H3  respectively  and  the  C2H2 
combines  with  the  oxygen  atom  to  produce  carbon  monoxide  (CO). 
As  a  whole,  under  very  low  oxygen  concentration,  in  which  all  of 
the  reaction  rates  weaken,  the  pathways  that  are  unimportant  at 
ordinary  oxygen  levels  now  could  become  prominent  and  grow 
in  relative  importance.  The  relevance  of  C2  channel  in  the  oxida¬ 
tion  kinetic  of  methane  under  MILD  condition  confirms  the  data  re¬ 
ported  in  the  results  of  some  Refs.  [15,20].  In  other  word 
importance  of  both  Low-Temperature  and  High-Temperature  reac¬ 
tion  pathways  in  oxidation  of  methane  under  MILD  condition 
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shows  the  kinetic  evolution  of  fuels  oxidation  process  is  different 
in  MILD  regime  from  traditional  flames. 

5.  Conclusion 

In  the  present  work,  methane  oxidation  to  CO  and  C02  is  stud¬ 
ied  under  MILD  condition  using  a  CFD  modeling  and  also  zero 
dimensional  analysis.  In  the  CFD  approach,  an  experimental  co¬ 
flow  burner  is  modeled  using  the  RANS  equations.  The  zero  dimen¬ 
sional  analysis  is  performed  for  a  WSR.  Results  show  that  under 
MILD  condition  chemical  pathways  including  conversion  of  methyl 
(CH3)  to  higher  hydrocarbons  (i.e.  ethane)  are  activated  and  some 
portions  of  methane  oxides  to  CO  through  conversion  to  heavier 
ethane  molecules.  That  means  that  methane  Low-Temperature 
reaction  pathways  activates  under  MILD  condition.  Moreover,  a  de¬ 
crease  of  02  concentration  in  oxidizer  does  not  always  leads  to 
lower  production  of  CO  and  C02.  This  can  be  because  of  dissocia¬ 
tion  of  C02  at  temperature  of  higher  than  1500  K.  Furthermore, 
the  ratio  of  rates  of  CO  producer  reactions  (R12  and  R13:  O  +  C2H2 
<^>  CH2(S)  +  CO)  in  ethane  oxidation  route,  to  CO  consumer  reaction 
(R59:  OFI  +  CO  <^>  H  +  C02)  in  the  main  route,  changes  with  chang¬ 
ing  of  02  level.  CFD  results  show  that  reaction  rates  of  R12  and  R13 
increase  and  become  larger  than  the  rate  of  R59  for  higher  02  con¬ 
centration  and  it  reveres  for  lower  02  concentration.  Also,  WSR 
analysis  shows  that  for  temperatures  lower  than  1500  K,  an  in¬ 
crease  in  temperature  increases  C02  concentration  and  decreases 
CO  concentration.  While  the  effect  of  temperature  on  CO  and  C02 
concentrations  reverses  for  temperature  higher  than  1500  K. 
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